Towards a global classification of excitable reaction— diffusion systems 
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Patterns in reaction-diffusion systems near primary bifur- 
cations can be studied locally and classified by means of am- 
plitude equations. This is not possible for excitable reaction- 
diffusion systems. In this Letter we propose a global classifi- 
cation of two variable excitable reaction-diffusion systems. In 
particular, we claim that the topology of the underlying two- 
dimensional homogeneous dynamics can be used to organize 
the system's behavior. 
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Many dissimilar experimental and model reaction- 
diffusion systems display similar behavior. This raises 
the question of whether a classification can be found so 
that seemingly unrelated systems with similar dynamics 
can be understood as belonging to the same equivalence 
class. We argue in this Letter that a global classification 
is possible for reaction-diffusion systems whenever the 
underlying homogeneous dynamics can be mapped onto 
a planar flow (i.e., a flow in R 2 ). This classification is 
done in terms of model families jl| that are largely de- 
termined by the homogeneous dynamics. In support of 
our conjecture we analyze experiments done in an open 
reactor using the FIS (Ferrocyanide-Iodate-Sulfite) re- 
action |^,D and two reaction-diffusion models that 
display similar patterns to those of the experiment, even 
though they are not accurate models of the FIS kinetics. 
We explain these common behaviors by noting that they 
all have similar homogeneous dynamics and discuss the 
main features of their model family. 

The classification of system behaviors lies at the heart 
of dynamical systems theory and here the normal form 
theorem is one of the main achievements |J. It states 
that close to a bifurcation point the dynamics of any 
sufficiently smooth system can be reduced to a simplified 
set of equations that is locally topologically equivalent to 
the full vector field [Q . The normal form approach is local 
in both phase and parameter space. Thus, classifications 
based on it provide generic descriptions only when the 
system is near threshold. For this reason, the patterns 
that occur in the excitable systems that we discuss in this 
Letter cannot be classified with normal form techniques. 

Excitability is a common dynamical behavior that oc- 
curs, for example, in the kinetics of neuronal membrane 



potentials M, in semiconductor lasers with feedback |9[| 
and in a variety of chemical reactions The picture 
of excitability we use is defined in terms of the spatially 
uniform dynamics. It assumes the existence of a stable 
fixed point (a spatially uniform stationary solution) such 
that perturbations above a threshold result in a large ex- 
cursion in phase space before the systems decays back 
to the fixed point. A spatially uniform excitable dynam- 
ics can result in wave propagation or in multistability 
when spatial variations are allowed. In particular this 
is true when the coupling is diffusive, i.e., for reaction- 
diffusion systems (see e.g. Since these behaviors 
arise from a finite perturbation of the spatially uniform 
steady state, classifications based on local normal forms 
cannot be used. 

We now consider a reaction-diffusion system whose ho- 
mogeneous dynamics can be mapped onto a planar flow. 
Because of the limited set of asymptotic behaviors al- 
lowed by flows in R 2 [[[o| they are particularly easy to 
classify. Given the fixed points of the flows there is a 
finite set of possible asymptotic behaviors and transi- 
tions among them JHJ]. Thus, we can classify families 
of planar flows that smoothly deform into one another 
by varying parameters |l2fl . It is important to note that 
this is a global description, as opposed to the local one 
provided by near threshold normal forms. This global 
description is valid for the spatially uniform dynamics. 
When diffusion terms are added the planar system be- 
comes high-dimensional. We conjecture that there ex- 
ists a natural global description for the high-dimensional 
reaction-diffusion system that is the direct analogue of 
the known global description of the two-dimensional pla- 
nar system. This means there should exist a model family 
of two coupled reaction-diffusion equations which dis- 
plays the same patterns and transitions among them as 
any other member of the class. The central statement 
is that, for reaction-diffusion systems which in the ho- 
mogeneous limit can be described by planar flows, these 
model families are mainly determined by the homogenous 
dynamics. 

Given a set of two reaction-diffusion equations with 
bounded diffusion terms, the local dynamics at each time 
and point in space is, properly understood, a perturba- 
tion of the homogeneous planar flow. Thus, the model 
family of the class which a particular system belongs to 
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should not only display homogeneous behaviors that are 
topologically equivalent to those of the system we are 
trying to classify, but it should also contain its closest 
perturbations. Since flows in R 2 with a finite number of 
fixed points have a finite set of asymptotic behaviors and 
can undergo a countable set of bifurcations, we then ex- 
pect to find a countable set of model families. This gives 
the desired classification of excitable systems describable 
by pairs of reaction-diffusion equations. In this Letter 
we present a particular application which supports our 
viewpoint. 




FIG. 1. Planar flow with one excitable fixed point (a) and 
its closest perturbations (b)-(d). 

We first develop a simple geometric model of excitabil- 
ity in R 2 suitable for our purposes. We consider a set 
of equations u — f(u,v), v — g(u,v), with only one at- 
tractor: a stable fixed point, x = (u,v). Excitability 
implies that the flow lines turn around before coming 
back to x. The simplest situation with this property is 
depicted in Fig. |l| (a), that corresponds to the case in 
which x (the square) is the only limit set of the flow. 
This flow can generically undergo a "saddle-repellor" , a 
saddle-node or a Hopf bifurcation, leading to the flows 
shown in Figs, [j] (b)-(d), respectively. The fixed point 
x is stable in all cases but in (d), where there is an at- 
tracting limit cycle. In Fig. |l| (b) and (c) it coexists with 
a saddle (the diamond) and a repellor (the white circle) 
or a node (the black circle). While the Hopf bifurca- 
tion requires that x be a spiral, which is unrelated to 
it being excitable, the saddle-node and saddle-repellor 
bifurcations may be linked to the excitability of x. In 
fact, excitability is related to the existence of a separa- 
trix, which becomes the stable manifold of the saddle 
after it is born. Also, the turn around of the flow lines 
implies that the orientation of the nullclines (the curves 
that satisfy f(u, v) = = g(u, v)) is such that, by not too 
large a perturbation, they can eventually become tangent 
at a point. Generically, when this happens, a saddle- 
node or saddle-repellor bifurcation occurs. These two 
types of bifurcations "meet" at a codimension two point, 
the Takens-Bogdanov point. In fact, the simplest fam- 
ily that contains the flows in Figs. [j] (a)-(c) outside a 



neighborhood of x is given by the Takens-Bogdanov nor- 
mal form u = v, v = fX\ + fj,2V + u 2 + uv. It is 
easy to modify these equations so that they also describe 
the behavior at and near x (see e.g. |J). Flows of this 
"extended" family with only one fixed point have the in- 
flection of the flow lines necessary for excitability. We 
may thus call it a general model of excitability in R 2 , for 
systems with one fixed point. Since the family contains 
all the relevant perturbations of the flow in Fig. [I] (a), it 
is also suitable in the extended case 
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FIG. 2. Transition from replicating spots to labyrinthine 
patterns in the FIS reaction when the homogeneous system 
approaches the saddle-node saddle-repellor bifurcation. 

We now consider a reaction-diffusion system with an 
underlying planar homogeneous dynamics that has only 
one stable but excitable fixed point. As a model family 
of its class we choose a two-variable reaction-diffusion 
system that, in the homogeneous limit, contains the 
flows in Figs, [j] (a)-(c). Now, any affine transformation 
(u, v) — y M.[u — ito, v — vq), with A4 a constant 2x2 ma- 
trix, will give another family with an equivalent homoge- 
neous dynamics and, in many cases, the right "winding" 
of the flow lines. If we introduce such a transformation 
in the set d t u = D u V 2 u + f(u, v), d t v — D v \7 2 v + g(u, v), 
we get a new set with cross-diffusion terms. Thus, once 
we have any family, u — f(u,v), v — g(u,v), with the 
"right" homogeneous dynamics, we expect, in most cases, 
that each (stationary) pattern of the system of interest 
be equivalent to a pattern of a flow in the family 

d t u = D uu V 2 u + D uv V 2 v + f(u, v), 

d t v = D vu V 2 u + D vv V 2 v + g(u, v). (1) 

Since it is possible to get rid of cross-diffusion by a simple 
transformation, it should always be possible to choose a 
particular model family without it. 

Now we discuss three concrete examples of systems 
with the excitable homogeneous dynamics discussed 
above. Consider the FIS reaction, which produces the 
patterns of Fig. ||. In the well-mixed (i.e., spa- 
tially homogeneous) case, it is known to exhibit excitabil- 
ity, bistability and oscillations |]lq| . These homogeneous 
behaviors can be described by planar dynamical sys- 
tems. The dynamical models of the system (a set of 10 
ODE's) |li| show that for the parameter values used in 
all the experiments, there is a separation of timescales 
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FIG. 3. A partial bifurcation set for the homogeneous 
Gray-Scott model (a) and two examples of patterns that oc- 
cur for the parameters indicated in (a) when D u = 1 and 
D v — 0.5. The line corresponds to the curve of saddle-repellor 
bifurcations. 

that allows the reduction of the original 10- 
dimensional system to a planar one. This planar sys- 
tem can have at most three fixed points, one of them a 
saddle. Experimentally, only stable solutions can be ob- 
served. The existence of the unstable fixed points and 
other limit sets are deduced from the model. 

In this Letter we discuss the replicating spots 
(Fig. |^ (a)) and lamellar structures (Fig. |^ (b)) that are 
found when there is only one homogeneous stationary 
solution, the low pH fixed point, and it is stable j|. Ex- 
perimentally, the spots are initiated by finite perturba- 
tions of the low pH state. This reflects its excitability. 
Thus, a model family like the one sketched before should 
reproduce at least some of the observed patterns. The 
transition from spots to lamellae shown in Fig. ^ is ob- 
served as the concentration of ferrocyanide, [Fe(CN) 6 ~], 
is decreased ||. The homogeneous system approaches a 
saddle-node bifurcation as [Fe(CN)g~] is decreased. Be- 
low a critical value there are three fixed points (a low 
pH one, a high pH one, and an intermediate pH saddle). 
Above this value only the low pH fixed point persists. 
We will show that both the patterns and this transition 
are contained in the model family we propose. 

The Gray-Scott and Fitzhugh-Nagumo models display 
spot replication and lamellar patterns. The Gray-Scott 
model is given by Q , jlTj : 

du 



.„ D u V 2 u-uv 2 + A(l-u), 
at 

dv n 9 

— = D V W v + uv z - Bv. 

dt 



(2) 



and it has been shown to behave similarly to the experi- 
ment g. The Fitzhgugh-Nagumo model Q, [§, can be 
written as: 



du 2 

— = D U V u- a(v + a x u - a ), 

— = D v V z v + v - v A + u. 
dt 



(3) 



The homogeneous dynamics of both models is a flow in 
R 2 . Furthermore, the families described by Eqs. (|J) and 
(||) contain the flows of Figs. [I] (a)-(c). Thus, both sys- 
tems have the "right" homogeneous dynamics to 




FIG. 4. Similar to Fig. |3[ but for the Fitzhgugh-Nagumo 
model. The examples in (b) and (c) were obtained for D u = 1 
and D v = 0.166. 

construct model families of one another and of the FIS 
reaction. We conclude that all these systems belong to 
the same equivalence class. Furthermore, for the qualita- 
tive comparison we present, cross diffusion terms are not 
necessary and either Eqs. (|J) or (|^) can be used as the 
model family for this class. 

We show in Figs. ||| and || snapshots of some of the pat- 
terns obtained in numerical simulations of Eqs. (|^) and 
(^), respectively, when the systems have only one stable 
but excitable fixed point. We observe spot replication 
and labyrinthine patterns. As in the experimental sys- 
tem, we observe a transition from spots to lamellar pat- 
terns as the homogeneous systems approach the saddle- 
node (saddle-repellor) bifurcation (see the partial bifur- 
cation sets shown in Figs. || (a) and [| (a)) pl|. We can 
provide a rough explanation for this if we think of sta- 
tionary solutions and regard diffusion as a perturbation 
on the homogeneous dynamics at each spatial point. In 
this sense, diffusion makes the system cross the saddle- 
node (or saddle-repellor) bifurcation. It is clear that the 
closer the homogeneous system is to the bifurcation, the 
smaller must be the perturbation and thus V 2 u and V 2 i>. 
For this reason, more spatially extended patterns with 
smaller Laplacians (such as lamellae) can be supported 
when the homogeneous system is closer to the bifurcation 
point. 

We have proposed a global classification scheme for 
excitable reaction-diffusion systems in terms of model 
families whose choice is based on their underlying homo- 
geneous dynamics. We have concluded that the FIS reac- 
tion and Eqs.(||) and (^), in the region of parameter space 
discussed, belong to the same equivalence class. This 
class is organized around a Takens-Bogdanov point ||(]] 
and we call it the Takens-Bogdanov model of excitabil- 
ity. Families with other types of homogeneous dynamics 
will serve as templates for other classes of excitable sys- 
tems. We believe that this classification approach will 
be possible whenever the homogeneous dynamics can be 
mapped onto a planar flow. The fact that limit sets of 
planar flows are so restricted may be the explanation 
of the ubiquitous presence of certain patterns in diverse 
systems. For example, this might be the reason that the 
complex Ginzburg-Landau equation reproduces behav- 
iors outside its range of applicability as a near-threshold 
normal form. 
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